Discovery of the first giant double radio relic in a galaxy cluster found in the PLANCK Sunyaev-Zel'dovich cluster survey: PLCK G287.0+32.9

We report the discovery of large scale diffuse non-thermal radio emission in PLCK G287.0+32.9, an exceptionally hot (T ~ 13 keV), massive and luminous galaxy cluster, strongly detected by the PLANCK satellite in a recent, all-sky blind search for new clusters through Sunyaev-Zel'dovich effect. Giant Metrewave Radio telescope 150 MHz and Very Large Array 1.4 GHz radio data reveal a pair of giant (>1 Mpc) "arc" shaped peripheral radio-relics (signatures of shock waves) of unprecedented scale (linear separation ~4.4 Mpc at redshift 0.39), located at distances from the cluster center which are about 0.7 and 1.3 of the cluster's virial radius. Another possible giant relic and a radio-halo is detected closer to the cluster center. These relic sources are unique "signposts" of extremely energetic mergers and shocks (both ongoing and past), that are assembling and heating up this very massive galaxy cluster. They are also a probe of the filamentary cosmic-web structure beyond the cluster virial radius. Optical imaging with the IUCAA 2 meter telescope and XMM-Newton X-ray data confirm a very rich galaxy cluster with a morphologically disturbed core region, suggesting a dynamically perturbed merging system.Comment: This version matches with the published version in Astrophysical Journal Letter

A 3D model for carbon monoxide molecular line emission as a potential cosmic microwave background polarization contaminant

We present a model for simulating Carbon Monoxide (CO) rotational line emission in molecular clouds, taking account of their 3D spatial distribution in galaxies with different geometrical properties. The model implemented is based on recent results in the literature and has been designed for performing Monte-Carlo simulations of this emission. We compare the simulations produced with this model and calibrate them, both on the map level and on the power spectrum level, using the second release of data from the Planck satellite for the Galactic plane, where the signal-to-noise ratio is highest. We use the calibrated model to extrapolate the CO power spectrum at low Galactic latitudes where no high sensitivity observations are available yet. We then forecast the level of unresolved polarized emission from CO molecular clouds which could contaminate the power spectrum of Cosmic Microwave Background (CMB) polarization B-modes away from the Galactic plane. Assuming realistic levels of the polarization fraction, we show that the level of contamination is equivalent to a cosmological signal with r 720.02. The Monte-Carlo MOlecular Line Emission (MCMole3D) Python package, which implements this model, is being made publicly available

MILCA, a Modified Internal Linear Combination Algorithm to extract astrophysical emissions from multi-frequency sky maps

The analysis of current Cosmic Microwave Background (CMB) experiments is based on the interpretation of multi-frequency sky maps in terms of different astrophysical components and it requires specifically tailored component separation algorithms. In this context, Internal Linear Combination (ILC) methods have been extensively used to extract the CMB emission from the WMAP multi-frequency data. We present here a Modified Internal Linear Component Algorithm (MILCA) that generalizes the ILC approach to the case of multiple astrophysical components for which the electromagnetic spectrum is known. In addition MILCA corrects for the intrinsic noise bias in the standard ILC approach and extends it to an hybrid space-frequency representation of the data. It also allows us to use external templates to minimize the contribution of extra components but still using only a linear combination of the input data. We apply MILCA to simulations of the Planck satellite data at the HFI frequency bands. We explore the possibility of reconstructing the Galactic molecular CO emission on the Planck maps as well as the thermal Sunyaev-Zeldovich effect. We conclude that MILCA is able to accurately estimate those emissions and it has been successfully used for this purpose within the Planck collaboration.Comment: 13 page

Far-Infrared to Millimeter Astrophysical Dust Emission. II: Comparison of the Two-Level Systems (TLS) model with Astronomical Data

In a previous paper we proposed a new model for the emission by amorphous astronomical dust grains, based on solid-state physics. The model uses a description of the Disordered Charge Distribution (DCD) combined with the presence of Two-Level Systems (TLS) defects in the amorphous solid composing the grains. The goal of this paper is to confront this new model to astronomical observations of different Galactic environments in the FIR/submm, in order to derive a set of canonical model parameters to be used as a Galactic reference to be compared to in future Galactic and extragalactic studies. We confront the TLS model with existing astronomical data. We consider the average emission spectrum at high latitudes in our Galaxy as measured with FIRAS and WMAP, as well as the emission from Galactic compact sources observed with Archeops, for which an inverse relationship between the dust temperature and the emissivity spectral index has been evidenced. We show that, unlike models previously proposed which often invoke two dust components at different temperatures, the TLS model successfully reproduces both the shape of the Galactic SED and its evolution with temperature as observed in the Archeops data. The best TLS model parameters indicate a charge coherence length of \simeq 13 nm and other model parameters in broad agreement with expectations from laboratory studies of dust analogs. We conclude that the millimeter excess emission, which is often attributed to the presence of very cold dust in the diffuse ISM, is likely caused solely by TLS emission in disordered amorphous dust grains. We discuss the implications of the new model, in terms of mass determinations from millimeter continuum observations and the expected variations of the emissivity spectral index with wavelength and dust temperature. The implications for the analysis of the Herschel and Planck satellite data are discussed.Comment: Accepted for publication in A&A (16 pages, 9 figures, 6 tables

Bubble size statistics during reionization from 21-cm tomography

The upcoming SKA1-Low radio interferometer will be sensitive enough to produce tomographic imaging data of the redshifted 21-cm signal from the Epoch of Reionization. Due to the non-Gaussian distribution of the signal, a power spectrum analysis alone will not provide a complete description of its properties. Here, we consider an additional metric which could be derived from tomographic imaging data, namely the bubble size distribution of ionized regions. We study three methods that have previously been used to characterize bubble size distributions in simulation data for the hydrogen ionization fraction – the spherical-average (SPA), mean-free-path (MFP) and friends-of-friends (FOF) methods – and apply them to simulated 21-cm data cubes. Our simulated data cubes have the (sensitivity-dictated) resolution expected for the SKA1-Low reionization experiment and we study the impact of both the light-cone (LC) and redshift space distortion (RSD) effects. To identify ionized regions in the 21-cm data we introduce a new, self-adjusting thresholding approach based on the K-Means algorithm. We find that the fraction of ionized cells identified in this way consistently falls below the mean volume-averaged ionized fraction. From a comparison of the three bubble size methods, we conclude that all three methods are useful, but that the MFP method performs best in terms of tracking the progress of reionization and separating different reionization scenarios. The LC effect is found to affect data spanning more than about 10 MHz in frequency (Δz ∼ 0.5). We find that RSDs only marginally affect the bubble size distributions

The ROSAT-ESO Flux Limited X-ray Galaxy Cluster Survey (REFLEX II) I. Newly identified X-ray luminous clusters at z>=0.2

We report 19 intermediate redshift clusters newly detected in the ROSAT All-Sky survey that are spectroscopically confirmed. They form a part of 911 objects in the REFLEX II cluster catalogue with a limiting flux of 1.8\times10^12 erg/s/cm2 in the 0.1-2.4 keV ROSAT band at redshift z >= 0.2. In addition we report three clusters from the REFLEX III supplementary catalogue, which contains objects below the REFLEX II flux limit but satisfies the redshift constraint above. These clusters are spectroscopically followed-up by our ESO NTT-EFOSC2 campaigns for the redshift measurement. We describe our observing and data reduction methods. We show how X-ray properties such as spectral hardness ratio and source extent can be used as important diagnostics in selecting galaxy cluster candidates. Physical properties of the clusters are subsequently calculated from the X-ray observations. This sample contains the high mass and intermediate-redshift galaxy clusters for astrophysical and cosmological applications.Comment: Astronomy and Astrophysics (in press

Non-thermal pressure in the outskirts of Abell 2142

Clumping and turbulence are expected to affect the matter accreted on to the outskirts of galaxy clusters. To determine their impact on the thermodynamic properties of Abell 2142, we perform an analysis of the X-ray temperature data from XMM-Newton via our SuperModel, a state-of-the-art tool for investigating the astrophysics of the intracluster medium already tested on many individual clusters (since Cavaliere, Lapi & Fusco-Femiano 2009). Using the gas density profile corrected for clumpiness derived by Tchernin et al. (2016), we find evidence for the presence of a non-thermal pressure component required to sustain gravity in the cluster outskirts of Abell 2142, that amounts to about 30 per cent of the total pressure at the virial radius. The presence of the non-thermal component implies the gas fraction to be consistent with the universal value at the virial radius and the electron thermal pressure profile to be in good agreement with that inferred from the SZ data. Our results indicate that the presence of gas clumping and of a non-thermal pressure component are both necessary to recover the observed physical properties in the cluster outskirts. Moreover, we stress that an alternative method often exploited in the literature (included Abell 2142) to determine the temperature profile kBT = Pe/nebasing on a combination of the Sunyaev-Zel'dovich (SZ) pressure Peand of the X-ray electron density nedoes not allow us to highlight the presence of non-thermal pressure support in the cluster outskirts

Planck 2013 results. XXII. Constraints on inflation

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to be ns = 0:9603 _ 0:0073, ruling out exact scale invariance at over 5_: Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0:11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V00 < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n _ 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining dns=dln k = 0:0134 _ 0:0090. We verify these conclusions through a numerical analysis, which makes no slowroll approximation, and carry out a Bayesian parameter estimation and model-selection analysis for a number of inflationary models including monomial, natural, and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for the post-inflationary entropy generation epoch, thus obtaining nontrivial data-driven constraints. We also present a direct reconstruction of the observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with second-order slow-roll predictions. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit by __2 e_ _ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with Planck bounds on fNL. Planck constrains with unprecedented accuracy the amplitude and possible correlation (with the adiabatic mode) of non-decaying isocurvature fluctuations. The fractional primordial contributions of cold dark matter (CDM) isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component can improve the _2 e_ by approximately 4 as a result of slightly lowering the theoretical prediction for the ` <_ 40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions

Impact on the tensor-to-scalar ratio of incorrect Galactic foreground modelling

A key goal of many Cosmic Microwave Background experiments is the detection of gravitational waves, through their B-mode polarization signal at large scales. To extract such a signal requires modelling contamination from the Galaxy. Using the Planck experiment as an example, we investigate the impact of incorrectly modelling foregrounds on estimates of the polarized CMB, quantified by the bias in tensor-to-scalar ratio r, and optical depth tau. We use a Bayesian parameter estimation method to estimate the CMB, synchrotron, and thermal dust components from simulated observations spanning 30-353 GHz, starting from a model that fits the simulated data, returning r<0.03 at 95% confidence for an r=0 model, and r=0.09+-0.03 for an r=0.1 model. We then introduce a set of mismatches between the simulated data and assumed model. Including a curvature of the synchrotron spectral index with frequency, but assuming a power-law model, can bias r high by ~1-sigma (delta r ~ 0.03). A similar bias is seen for thermal dust with a modified black-body frequency dependence, incorrectly modelled as a power-law. If too much freedom is allowed in the model, for example fitting for spectral indices in 3 degree pixels over the sky with physically reasonable priors, we find r can be biased up to ~3-sigma high by effectively setting the indices to the wrong values. Increasing the signal-to-noise ratio by reducing parameters, or adding additional foreground data, reduces the bias. We also find that neglecting a 1% polarized free-free or spinning dust component has a negligible effect on r. These tests highlight the importance of modelling the foregrounds in a way that allows for sufficient complexity, while minimizing the number of free parameters.Comment: 11 pages, 7 figures, submitted to MNRA